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Abstract. The process of energy release in the magnetosphere of a fast rotating, magnetized white dwarf can be explained in 
terms of the canonical spin-powered pulsar model. Applying this model to the white dwarf companion of the low mass close 
binary AE Aquarii leads us to the following conclusions. First, the system acts as an accelerator of charged particles whose 
energy is limited to £ p S= 3 TeV and which are ejected from the magnetosphere of the primary with the rate L^ m 5= 10 32 erg s -1 . 
Due to the curvature radiation of the accelerated primary electrons the system should appear as a source of soft y-rays (~ 
100 keV) with the luminosity < 3 x 10 27 erg s . The TeV emission of the system is dominated by the inverse Compton scattering 
of optical photons on the ultrarelativistic electrons. The optical photons are mainly contributed by the normal companion and 
the stream of material flowing through the magnetosphere of the white dwarf. The luminosity of the TeV source depends on the 
state of the system (flaring/quiet) and is limited to < 5 x 10 29 erg s . These results allow us to understand a lack of success in 
searching for the high-energy emission of AE Aqr with the Compton Gamma-ray Observatory and the Whipple Observatory. 
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1. Introduction 



As shown by Usov ( 1988 1, the process of energy release in the 
magnetosphere of a fast rotating, magnetized white dwarf can 
be explained in terms of a canonical spin-powered pulsar model 
(e.g. Arons & Scharlemann 1979) provided its surface temper- 
ature is r(i? w( j) % 10 6 K. The appearance of such a white dwarf 
is similar to that of a spin-powered neutron star (i.e. a canonical 
spin-powered pulsar) in several important aspects. In particular, 
its energy budget is dominated by the spin-down power 



L sd - 3 x 10 35 /4 



(—)' 
UOs/ 



ergs 



(1) 



which is spent mainly for the ejection of a relativistic wind. 
Here P s is the spin period of the white dwarf and fj.34 is its 
dipole magnetic moment expressed in units of 10 34 G cm 3 . The 
pressure of the wind significantly exceeds the ram pressure of 
the material surrounding the star at the accretion (Bondi) ra- 
dius and the particles accelerated in its magnetosphere reach 
the TeV energies (Usov 1988). According to Lipunov dl992> . 
the state of the white dwarf can be classified as an ejector. 
Following this we will refer such objects as Ejecting White 
Dwarfs (EWDs). 
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The relativistic particles accelerated in the magnetospheres 
of EWDs manifest themselves in the high-energy part of the 
spectrum (mainly due to the curvature radiation and the inverse 
Compton scattering of thermal photons), while the dissipation 
of the back-flowing current, which closes the current circuit in 
the magnetosphere, leads to a local heating of the white dwarf 
surface in the magnetic pole regions. Therefore, EWDs are ex- 
pected to appear as non-thermal y-ray pulsars and pulsing X- 
ray/UV sources with a predominantly thermal spectrum (see 
Usov lT995l . 

A detailed analysis of the origin and appearance of EWDs 
has been performed in the frame of the interpretation of non- 
identified galactic y-ray sources by Usov (1988), and of the 
anomalous X-ray pulsar IE 2259+586 by Paczyhski ( 1990) and 
Usov J1993I IT994I . As they have shown, the basic appearance 
of these objects fits well into the EWD-model and the forma- 
tion of fast rotating, magnetized white dwarfs should be a rel- 
atively common phenomenon in our Galaxy. Nevertheless, fur- 
ther development of these approaches was not effective mainly 
because of a lack of evidence for the white dwarf nature of 
these sources (see e.g. Mereghetti et al. [2002 and references 
therein). 

A target, to which the application of the EWD-model is cur- 
rently under consideration, is the degenerate component of the 
low mass close binary system AE Aqr. This star is unambigu- 
ously identified with a fast rotating (P s ^ 33 s) white dwarf 
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(Eracleous et al. 1994 > whose spin-down power exceeds its lu- 
minosity by at least a factor of 120 (de Jager et al. 1994 , Welsh 
1999 1. The appearance of the system significantly differs from 
that of magnetic cataclysmic variables whose radiation is pow- 
ered by the accretion process (Warner 1995). Instead, AE Aqr 
ejects material in the form of a non-relativistic stream (Welsh 
et al. 119981 Ikhsanov et al. [2004 1, and relativistic particles re- 
sponsible for the radio and, possibly, high-energy emission of 
the system (de Jager ll99"4V 

In this paper we analyze the appearance of AE Aqr in the 
high energy parts of the spectrum. The basic information on 
the system and, in particular, its appearance in TeV y-rays are 
briefly outlined in Sect. 2. The process of particle acceleration 
in the magnetosphere of the white dwarf within the EWD- 
model is discussed in Sect. 3. In Sections 4 and 5 we address 
the appearance of the leptonic and hadronic components of the 
accelerated particles, respectively. A brief summary and dis- 
cussion of acceleration mechanisms alternative to the EWD- 
model are given in Sect. 6. 



2. AE Aquarii 

AE Aqr is a non-eclipsing binary system with an orbital period 
of 9.88 hours and eccentricity of 0.02, which is situated at the 
distance of 100 pc. The system inclination angle and the mass 
ratio are limited to 48° < i < 62°, and 0.58 5 (q = M 2 /M{) < 
0.89 (Welsh et al. [T995l . 

The system emits detectable radiation in almost all parts 
of the spectrum. It is a powerful non-thermal radio source 
(Bastian et al. 1988 1 and, possibly, a y-ray emitter (see below). 
Its optical, UV, and X-ray radiation is predominantly thermal 
and comes from at least three different sites. The visual light 
is dominated by the normal companion (secondary), which is a 
K3-K5 red dwarf on or close to the main sequence (Welsh et al. 
fT995l .The primary is a fast rotating, magnetized white dwarf. 
The remaining light comes from a highly variable extended 
source, which manifests itself in the blue/UV continuum, the 
optical/UV broad single-peaked emission lines, and the non- 
pulsing X-ray component. This source is associated with the 
stream of material, which flows into the Roche lobe of the white 
dwarf from the normal companion, interacts with the magne- 
tosphere of the primary and is ejected out from the system 
without forming a disk (Wynn et al. 119971 Welsh et al. 119981 



Ikhsanov et al. 2004). This source is also suspected of being 
responsible for the peculiar rapid flaring of the star (Eracleous 
& Horne [T996l . 

2.1. Parameters of the white dwarf 

A justification of the white dwarf nature of the primary has 



been presented by Eracleous et al. (1994) on the basis of 
HST observations. As they have shown, the optical/UV puls- 
ing emission comes from two hot spots with the projected area 
A p ~ (0.8 - 4.3) x 10 16 cm 2 and the temperature 20000 
T po i £ 47000 K. Associating these spots with the magnetic pole 
regions, they have limited the angle between the rotational and 



magnetic axes of the primary to 75 % [3 % 11 and evaluated 
the temperature of the rest of its surface as 10000-16000 K. 

The white dwarf is spinning-down with a mean rate Po = 
5 .64 x 1 0~ 14 s s~ 1 , which implies its spin-down power (de Jager 
et al. lT994l Welsh [T999t 



L sd = 6 x 10 33 7 50 



(33 s) [pj 



erg s 



(2) 



Here /50 is the moment of inertia of the white dwarf expressed 
in units of 10 50 g cm 2 . 

L S d exceeds the luminosity of the system observed in the 
UV and X-rays by a factor of 120 and its persistent bolometric 
luminosity by a factor of 5. As mentioned above, such a situa- 
tion is typical for EWDs. The dipole magnetic moment of the 
white dwarf within this approach can be evaluated as (Ikhsanov 
1998) 



(to 



1.4 x 10 J 



(— r 

\33s/ 



6 x 10 33 ergs - 



1/2 



G cm 



(3) 



This implies that the strength of the surface magnetic field of 
the white dwarf in the magnetic pole regions is (Alfven & 
Falthammar H963l 



Bo = 



3-3 



^- 100 Rio 

R 3 8.8 

WQ 



1.4 x 10 34 Gcm 3 



MG, 



(4) 



and, correspondingly, the surface field strength at its magnetic 
equator is Bq/2 = 50 MG. Here Z?g.g is the radius of the white 
dwarf expressed in units of 10 8 8 cm. 

As recently shown by Ikhsanov et al. (2004), the above 
mentioned estimates contradict none of the currently observed 
properties of the system, and the Doppler Ha tomogram simu- 
lated within this approach is in a very good agreement with the 
tomogram derived from spectroscopic observations of AE Aqr. 
On this basis estimates (|3} and (0J will be used in our analysis. 

2.2. High-energy emission 

A search for TeV emission of AE Aqr was performed by three 
independent groups. The Potchefstroom group has reported on 
about 3 10 hours of observations in 1988-93 (Meintjes et al. 
119921 11994> . A few events of pulsing TeV emission with an 
averag ed flux (0.2 - 2) x 10~ 10 cnr 2 s _1 (at a threshold energy 
2.4 TeV) were detected during this period. Additionally, they 
reported a detection of two short (1 and 3 minutes duration) 
unpulsed outbursts on 25 June 1993 (the orbital phases 0.04 
and 0.05). The flux of the TeV emission during these events 
was about 8 x 10~ 9 cirr 2 s~\ which, under the assumption of 
isotropic emission, corresponds to the source luminosity of 
1 Here d\oo is the distance to AE Aqr nor- 



10 34 d*,,, ergs" 



malized to lOOpc. 

A detection of TeV emission from AE Aqr has also been 
reported by the Durham group (Bowden et al. 1992 Chadwick 
et al. [T995l . They recorded about 170 hours of data in 1990-93. 
A persistent component of pulsing (at 60.46 mHz) TeV emis- 
sion with a time-averaged flux of about 10~ 10 cirT 2 s _! (at en- 
ergy > 350 Ge V) was detected. Additionally, they reported the 
detection of a 1 minute (13 October 1990, orbital phase 0.40) 
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and a 70minutes (11 October 1993, orbital phases 0.62-0.74) 
outburst. The luminosity of the y-ray source during these out- 
bursts was ~ 6 x 10 33 d 2 m) erg s" 1 and 3 x 10 32 ii 2 00 erg s" 1 , re- 
spectively. In both cases the signal was modulated with a half 
of the spin period of the white dwarf (60.46 mHz). 

On the other hand, the analysis of data recorded by the 
Whipple group (68 hours of observations in 1991-95) have 
shown no evidence for any steady, pulsed or episodic TeV 
emission of AE Aqr (Lang et al. 11998b . The upper limits to 
the flux of a steady emission, $ 4.0 x 10" 12 cm" 2 s" 1 , and a 
pulsing emission S 1.5 X 10" 12 cm" 2 s" 1 , at an energy threshold 
of 900 Ge V have been derived. This suggests that the persistent 
luminosity AE Aqr in TeV energy range is unlikely to exceed 
5 x 10 30 erg s . At the same time, the database of the Whipple 
group is too small for any conclusions on the non-frequent TeV 
y-ray events reported by the Potchefstroom and Durham groups 
(for a discussion see Lang et al. U998> . 

Finally, a negative result of a search for 0.1-1 GeV emis- 
sion from AE Aqr with the Compton Gamma-Ray Observatory 
(CGRO) has been reported by Schlegel et al. fH)9"5l The de- 
rived upper limits suggest that the luminosity of the system in 
the EGRET energy range is smaller than 10 31 erg s _1 . 

A lack of success in searching for y-ray emission from 
AE Aqr reported by the EGRET and Whipple groups raises 
the following questions: Do we have any theoretical grounds to 
suggest that this system can be an emitter of high-energy radi- 
ation and particles and if so, how bright this source could be ? 
In order to answer these questions we performed an analysis of 
particle acceleration based on the EWD-model of the system. 
We show that the expected intensity of high-energy emission 
within this model is indeed below the threshold of detectors 
used by the CGRO and Whipple observatories. 



3. Particle acceleration 

We consider a process of particle acceleration in the magneto- 
sphere of a white dwarf, whose rotation period P s = 33 P33 s 
(the angular velocity Q = 2n/P s =t 0.19 rads" 1 ), the mass 
M wd = 0.8 M .8 Mq, the strength of the magnetic field at the 
magnetic pole regions Bq = 10 8 G and the surface tempera- 
ture T < 50000 K. Under these conditions a scale height of its 
atmosphere, 

is significantly smaller than both the radius of the white dwarf 
(^wd - 10 8 ' 8 cm) and the lower limit to the radius of its polar 
caps, 

Afl p - l^j'Vd * 4 x 10 7 < 8 2 (|^) cm. (6) 

This implies that the vacuum approximation to the medium sur- 
rounding the white dwarf is applicable and therefore, the pro- 
cess of particle acceleration can be treated in terms of EWD- 
model. 

According to Arons & Scharlemann ( 1979 1, the component 
of the electric field E\\ = (E ■ B)/\B\ along the magnetic field 



B, which is generated in the polar cap regions of a fast rotating 
magnetized star surrounded by a vacuum can be evaluated as 



E« - E AS x 




at < s < M p , 
at AR p < s < R wd , 



V2/W" at s > R wA , 



where 



'-AS 



1 (£W 



wd 



8 V3 I 



5/2 



Bo, 



(7) 



(8) 



and r = (R w a + s) is a distance from the surface of the white 
dwarf. This means that the electric potential in the magneto- 
sphere of the white dwarf is 

tpasir) = I E\\ds - 2V2£ , AS^wd 



jr. 



r 

B w d 



1 



(9) 



The total electric potential in the magnetosphere of a spin- 
powered neutron star contains two additional terms which rep- 
resent the effects of general relativity (Muslimov & Tsygan 
1992), and the effect of e ± pairs creation in the outer magne- 



tosphere (Cheng et al. 1986). However, a contribution of these 
terms to the electric potential in the magnetosphere of a white 
dwarf (whose radius exceeds that of a neutron star by a factor 
of 500 and the magnetic field is too weak for a pair creation 
process to be effective) are small and can be neglected. 

An application of the EWD-model to the white dwarf in 
AE Aqr implies a validity of the following two assumptions. 
The first is a lack of accretion of material onto the surface 
of the white dwarf. If this assumption were not valid the star 
would appear as an accretion-powered X-ray pulsar. The sys- 
tem is indeed an emitter of X-rays. But all of currently estab- 
lished properties of the emission are inconsistent with those 
predicted by the accretion-powered pulsar model. In particu- 
lar, the X-ray spectrum of AE Aqr is soft and significantly dif- 
fers from the typical spectra of accretion-powered white dwarfs 
(Clayton & Osborne [T995l Choi et al. U999> . Furthermore, the 
surface temperature of the white dwarf at the magnetic pole re- 
gions evaluated by Eracleous et al. (1994) is < 50000 K, i.e. 
a factor of 2000 smaller than the typical temperature at the 
base of an accretion column. Finally, observations of AE Aqr 
with XMM-Newton recently reported by Itoh et al. (2005 1 sug- 
gest that the number density of plasma responsible for the de- 
tected X-rays is n e ~ 10 11 cm" 3 (i.e. a few orders of magnitude 
lower than corresponding conventional estimates in the post- 
shock accretion column) and that the linear scale of the source 
is i v J 2 x 10 I0 cm (i.e. a factor of 40 larger than the radius 
of the white dwarf). This clearly indicates that the source of 
the observed X-ray emission is located at a significant distance 
from the white dwarf and therefore, cannot be powered by an 
accretion of material onto its surface. 

The second assumption is that the material streaming 
through the magnetosphere of the white dwarf does not sig- 
nificantly effect the electric potential generated at its surface. 
A hint for a validity of this assumption comes from the model 



of Michel & Dessler (1981) who have considered a situation 
in which a spin-powered pulsar is surrounded by a dead disk. 
According to their results a presence of a disk inside the light 
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cylinder of the pulsar would not suppress the electric potential 
at the stellar surface. Instead, an interaction between the field 
lines and the disk material leads to an additional term in the ex- 
pression of the total electric potential at the stellar surface. This 
finding has already been used by de Jager ( 1994 > for construct- 
ing a dead disk model for AE Aqr. The value of the electric 
potential in this case is comparable to that of <p as (ro), where ro 
is the radius of the disk (or a distance of closest approach of 
the stream to the white dwarf). In this light expression (|9j rep- 
resents an upper limit to the electric potential which could be 
generated in the magnetosphere of the white dwarf within the 
EWD-model. 

The energy of particles accelerated in the (f as {r) potential 
can be limited to £ p ^ £ , where 



4. Radiative losses of primary electrons 

The radiative losses of ultrarelativistic electrons accelerated in 
the magnetosphere of a magnetized compact star are governed 
mainly by two mechanisms: i) the curvature radiation, and 
ii) the inverse Compton scattering of thermal photons (Michel 
19911. 



4.1. Curvature radiation 

The intensity and the mean photon energy of the curvature radi- 
ation emitted by ultrarelativistic electrons in the magnetosphere 
of the white dwarf are (Usov 1988) 

2e 2 c 



£ -2xl0 11 J U34.2 ^ 



( — ) 

\33s/ 



-5/2 



1/2 



1 



" 3 lT 



eV. (10) 



3 hey 



Here Iq is the scale of the acceleration region. A precise de- 
termination of this parameter is an open problem so far. In 
general case the value of Iq can be limited to (e.g. Beskin 
et al. fT993l (tf wd + AR p ) < l < R lc , where R ic = c/Q * 
1.6 x 10 11 (P/33 s)cm is the radius of the light cylinder of the 
white dwarf. Following this limitation one can estimate the en- 
ergy of the accelerated particles as 

8 GeV < £ p < 3 TeV. (11) 

The kinetic luminosity of the particle beam can be evalu- 
ated as 

iki„ * eip.Jh)N, (12) 
where 

N = 7r(Afl p ) 2 « GJ 0R wd )c (13) 
is the flux of ultrarelativistic particles from the white dwarf and 
(O • B) , ,_ A I P s r 1 / B () 



2 R c 



(19) 



(20) 



Here R c (Rc/Q.) 1 ^ 2 is the curvature radius of the magneto- 
spheric field lines and y = & p /mc 2 is the Lorentz factor of 
electrons which in the case of AE Aqr can be expressed using 
Eq. GDJ as 

p s v-5/2 \l . \l/2 



7e(r) 



4xiow 



- 



- 1 



(21) 



Since the Lorentz factor of the electrons and the values of the 
parameters / c (r) and E cmv (r) increase with r, the largest fraction 
of the curvature radiation is generated near the light cylinder of 
the white dwarf. Therefore, the upper limit to the mean photon 
energy and the intensity of the curvature radiation (L curv - l c N) 
expected from AE Aqr can be evaluated from Eqs. dl9l l - ( 12 li 
as 

y e (R lc ) * 6 X 10 6 



«Gj(^wd) 



2nce 



5 x 104 (sil) (wq) cm 3 (14) £curv(/?lc) * 100 keV 



is the Goldreich- Julian density above the primary surface. 

Combining Eqs. and ij^Ji with Eqs. fl!2i — J14i and setting 
k) = R\c yields 

1 Cl 4 R b A Bl 

L£n X = wd_J) cQ ^ Q 04Lsd) (15) 

4 V6 c 3 
where it is taken into account that 

L sd = f — =a_2. (16) 

3 c J 

Thus, the maximum kinetic luminosity of the beam of ultrarel- 
ativistic particles ejected by the white dwarf in AE Aqr is 

£sd 



icurv(^ic) * 3 x IOV34.2 (3^) er § s 



L--2X10 



32 



ergs 



(17) 



6x lO^ergs" 1 , 

It is comparable to the bolometric luminosity of the system dur- 
ing its optical quiescent state (see Beskrovnaya et al. U996t . 

If, however, the scale of the acceleration region is close to 
its minimum value, i.e. Iq ~ (R w <i + A/? p ), the kinetic luminosity 
of the beam of GeV-electrons does not exceed 



Hence, the curvature radiation of the primary electrons con- 
tributes to the system emission mainly in the soft (~ 100 keV) 
y-rays with a mean flux of ~ 10 _I5 ergcm _2 s _I . This value 
is significantly smaller than the threshold sensitivity of detec- 
tors used in the EGRET experiment. This provide us with a 
reasonable explanation of unsuccessful search for high-energy 
emission of the system with CGRO. 

4.2. Inverse Compton scattering 

The mechanism, which can be responsible for the very high- 
energy emission of the system, is the inverse Compton scat- 
tering of the thermal photons on the ultrarelativistic electrons. 
The energy of the scattered photons is as high as m e c 2 y(R\ c ) ~ 
3 TeV, where m e is the electron mass. 

The intensity of the TeV radiation generated due to the 



2y 



6 x 10 33 ergs- 



erg s 



(18) 



inverse Compton scattering is (Ochelkov & Usov 1983J the 
Klein-Nishina cross-section case) 



k s N, 



(22) 
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where 

3 / //?,.(- 



6fl, 



In 



46 th y 

OT e C 2 



11 



(23) 



is the rate of energy losses of an electron. Here e t h - 3kT is the 
mean energy of the thermal photons, <x T is the Thomson cross 
section, k is the Boltzmann constant, T is the mean temperature 
of the thermal emission and 



uq ^ al 



(Re 



(24) 



is the radiation energy density. Rq denotes the effective ra- 
dius of the source of the thermal emission, r s is the mean 
distance from this source to the region of scattering and a = 
7.6 X 10~ 15 erg cm -3 deg 4 is the radiation density constant. 

The field of the thermal radiation in the magnetosphere of 
the white dwarf in the case of AE Aqr is contributed by three 
separate sources: the normal companion, the white dwarf and 
the stream of material flowing through its magnetosphere. The 
relative contribution of these sources depends on the location 
of the scattering region and the state of the system (quiescent 
or flaring). As we require the energy of the scattered photons to 
be ~ TeV, the closest distance of the region of their generation 
to the surface of the white dwarf is r ^ 2x 10 10 cm (see Eq. l21> . 
The contribution of the white dwarf to the radiation field at this 
distance is significantly smaller than those of the normal com- 
ponent and of the stream. Furthermore, the ratio 77 = l cs /eE\f 
in the vicinity of the white dwarf is 



This indicates that the thermal emission from the surface of the 
white dwarf is not essential for the process of particle accel- 
eration. On this basis we neglect the contribution of the white 
dwarf in further consideration. 



As recently shown by Ikhsanov et al. (2004), the closest 
approach of the stream to the surface of the white dwarf within 
the EWD-model is limited to 



ro 



> 4 x 10 10 770.37 lxf A2 M" 1 / 7 M- 7 2/7 cm, 



(25) 



where M\-j is the mass transfer rate expressed in units of 
10 17 gs~' and 770.37 = rj/Q.31 is the parameter accounting for 
the geometry of the accretion flow, which in the case of a 
stream is normalized following Hameury et al. ( 1986). This in- 
dicates that the contribution of the stream to the radiation field 
during the quiescent state (M < 10 17 gs~') is comparable to 
that of the normal companion inside the magnetosphere and 
is almost negligible at the distances close to the light cylin- 
der (R\ c = 1.57 x 10 11 cm). Therefore, using the parameters of 
the secondary (i.e. T e s ^ 4500 K, Ro ^ Rq) one can evalu- 
ate the intensity of the TeV emission produced by the inverse 
Compton scattering during the quiescent state of the system as 
(seeEqs.|22H24} 

(26) 



I± <, 3 x 10 27 ergs 



Here we assumed that the scattering occurs at the light cylinder 
of the primary. 



The most favorable conditions for the generation of TeV 
emission due to the inverse Compton scattering occur during 
the flaring state. According to Beskrovnaya et al. {19961, the 
system emission during the strongest optical flares (which last 
a few minutes) is dominated by a source with an effective tem- 
perature ~ 20000 K and an effective area of ^ 10 20 cm 2 . As this 
source is situated inside the magnetosphere (at the distances 
of about ro) its contribution to the radiation field significantly 
exceeds that of the normal companion (the optical luminosity 
of the system during these events reaches ^ 10 33 ergs _1 ). In 
this case the generation of TeV emission due to the inverse 
Compton scattering reaches its highest efficiency in a region 
situated at the distances of about ro from the surface of the pri- 
mary and its intensity (according to Eqs.l22l- I24> can be limited 
to 



Z4 < 4 x 10 29 ergs _1 



(27) 



(the mass transfer rate during the strong flares according to the 
EWD-model is ~ 5 x 10 17 g s" 1 ). 

Thus, the upper limit to the intensity of TeV radiation of 
AE Aqr produced due to the inverse Compton scattering of 
thermal photons on the ultrarelativistic electrons is a factor 
of a few smaller the threshold sensitivity of detectors used in 
observations of the system with the Whipple y-ray telescope. 
On the other hand the evaluated luminosity of the source is 
high enough for its emission to be detected in the MAGIC, 
VERITAS and HESS experiments. 

5. Appearance of a proton beam 

One can also envisage a situation in which a significant frac- 
tion of L s d is converted to ultrarelativistic protons. The upper 
limit to the kinetic luminosity of the proton beam in this case 
is expressed by Eq. d!2i . The radiative losses of the protons 
in the magnetosphere of the white dwarf (due to the curvature 
radiation and the inverse Compton scattering) are significantly 
smaller than those of electrons and in the first approximation 
can be neglected. The energy of the beam, however, can be 
effectively converted into VHE y-rays if the trajectories of pro- 
tons intersect a target of a relatively dense (30 - 80 cirT 2 ) back- 
ground material. In this case a creation of 7r-mesons (n°, 7r ± ) 
and their subsequent decays into two VHE y-photons with the 
energy of about the energy of the primary relativistic protons 
would be expected. 

Two targets with the required column density can be indi- 
cated in AE Aqr. Namely, the atmosphere of the normal com- 
panion, and the stream of material moving through the mag- 
netosphere of the white dwarf. The interaction of TeV protons 
with these targets will produce a flux of TeV photons. However, 
can this emission be detected by an observer at the Earth ? A 
positive answer on this question in the light of currently estab- 
lished geometry of the system is not obvious. Indeed, the pho- 
tons generated in this process are strongly beamed along the 
incident particle's velocity vector (see e.g. Vestrand & Eichler 
1982 1. Therefore, an observer can detect only those photons 
that are produced when protons streaming toward him strike 
intervening target material. This means that the flux of the TeV 
photons generated in the atmosphere of the secondary can be 



6 



N.R. Ikhsanov and P.L. Biermann: High-energy emission of fast rotating white dwarfs 



detected only at the orbital phases (2n — A(p) ^ ip iS (2n + Aip), 
where 



Aip - — arctan I 
r 2n \ A 



1 



■arctan (0.378 ? - a2084 ), 



and if the inclination angle of the system is 
/^/ - - -arctan (0.378 ? - a2084 ). 



(28) 



(29) 



Here Rui) is the volume radius of the Roche lobe of the sec- 
ondary and A is the orbital separation. For the parameters of 
AE Aqr conditions {28} and {29) are 0.94 <, <p < 1 .06 and 

o o 

i > io ~ 67 - 69 , respectively. Hence the atmosphere of the 
secondary is definitely not a site of TeV emission detected at 
orbital phases 0.07 - 0.93. Furthermore, it is unlikely a site of 
the TeV emission at the rest of the orbital phases since the as- 

o 

sumption i > 67 contradicts optical observations of the system 
(see Sect.|2j. 

The TeV emission generated due to the interaction between 
the relativistic protons and the stream could be observed in a 
significantly wider range of orbital phases. However, the con- 
dition to the inclination angle of the system in this case is 



i > ; ~ arctan — 

2 \r 



(30) 



where r str is the radius of the stream cross section, which in 
the considered case is limited to (see e.g. Wynn et al. 119971 
Ikhsanov et al. l2004> : 



r str <(Q L1 /n) l/2 ^ 1.6x10' 



1 4500 K/ 



cm, 



(31) 



and ro is the distance of closest approach of the stream to the 
white dwarf, which in the case of EWD-model is expressed by 
Eq. d25i . Here Gli is the effective cross section of the mass 
transfer throat at the Lagrange point LI and T e ff is the effective 
temperature of the photosphere of the secondary. Combining 
Eqs. d25i . (13 0> and (13 1 i one finds that the TeV photons gener- 
ated due to interaction between the accelerated protons and the 
stream could be detected by an observer if the inclination angle 

o 

of the system were <: 86 . This, however, is inconsistent with 
the value of i derived from observations (otherwise the system 
would be an eclipsing binary). 

Finally, we would like to point out that in case of any targets 
a strong collimation of the beam of ultrarelativistic protons is 
required for the effective production of TeV emission. Indeed, 
the reported luminosity of the TeV source during outbursts is 
an order of magnitude larger than the upper limit to the ki- 
netic luminosity of the beam evaluated in Sect.[3] Taking into 
account that the efficiency of the energy transfer from protons 
to 7r°-mesons does not exceed 20% (Ozernoy et al. 1973 1, and 
that only a half of the energy of the 7r°-mesons is transferred to 
y-rays detected by an observer one can limit the opening angle 
of the beam to 



AT < 0.1% 



^kin 



9xl0 32 ergs-V\10 34 ergs 



rad, 



where 



x 



I 

ft/Gb 



for 
for 



Gb * Qt, 
Gb > ft- 



(32) 



(33) 



Here Gb and Q t are the cross sections of the beam and the target 
in the direction, perpendicular to the velocity vector of the rela- 
tivistic protons. This means that the opening angle of the beam 

o 

for any targets should be less than 7 . This value corresponds 
to the opening angle of the dipole magnetic field at the distance 
r\ - 2.5 7? w d from the surface of the primary. Therefore, to meet 
this criterium one has to assume that either a target is situated 
very close to the surface of the white dwarf or the accelerated 
protons at r > r\ do not follow the magnetic field lines but 
propagate through the magnetosphere along almost straight tra- 
jectories. Both of these assumptions, however, are inconsistent 
with either the predictions of the EWD-model or the currently 
reconstructed picture of the mass-transfer in AE Aqr. 

6. Summary and discussion 

Application of the EWD-model to AE Aqr shows that the high- 
energy emission of the system is dominated by the radiative 
losses of TeV electrons accelerated in the magnetosphere of 
the primary. The energy of these electrons is converted into 
the TeV photons mainly due to the inverse Compton scattering. 
The luminosity of the TeV source depends on the optical state 
of the system and lies within the interval (3-500)x 10 27 erg s _1 . 
Therefore, the expected maximum flux of TeV radiation from 
AE Aqr (during the strongest optical flares) within our model is 
limited to 55 2 x 10~ B cirT 2 s . This is below the threshold of 
detectors used in all TeV observations of the system reported 
so far and therefore, a lack of success in searching for the high- 
energy emission of AE Aqr by the Whipple group proves to be 
quite understandable. 

On the other hand, the origin of TeV emission detected by 
the Potchefstroom and Durham groups within our approach ap- 
pears to be a miracle. If this radiation is indeed emitted by 
AE Aqr one has to assume that an acceleration mechanism, 
which is more powerful than that investigated in this paper, op- 
erates in the system. What kind of a mechanism could it be ? 

A lack of a disk in the system forces us to reject a number 
of previously suggested acceleration scenarios. In particular, 
the unipolar inductor model of Cheng & Ruderman ( 1991 1 and 
the dead disk model of Michel & Dessler J198H are not appli- 
cable in the light of current view on the mass-transfer process 
in AE Aqr. 

The statistical acceleration mechanisms operating inside 
the magnetosphere of the primary are not effective due to a rel- 
atively small scale of the system. As shown by Kuijpers et al. 
( 1997 1, the typical Lorentz factor of electrons accelerating by 
the magnetic pumping in the magnetosphere of the white dwarf 
is about 100. This indicates that the statistical mechanisms of 
acceleration can be helpful for the interpretation of the radio 
emission of the system but their contribution to the very high 
energy y-rays is negligible. 

One can also envisage a situation in which the system acts 
as an unipolar inductor due to the interaction between the mag- 
netic field of the primary and the normal companion, which is 
partly situated inside the light cylinder of the white dwarf. This 
interaction leads to a perturbation of the magnetospheric field 
broadening the area of the hot spots at the magnetic pole re- 
gions by a factor of (^i c /^li) 1/2 - 1 -25 and creating an electric 
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potential in the region of interaction (here is the distance 
from the white dwarf to the LI point). The rate of energy re- 
lease due to this interaction is limited to 

Q2CB 2 (R L1 ) 33 2 6 / R 2 \ 2 _j 
2 frr " 10 ItTTo^) ei ' gS ' (34) 

where Qo is the effective cross section of the interaction be- 
tween the normal companion of the radius Ro and the magnetic 
field of the primary. As seen from Eq. (I34> . the total power of 
this interaction is at least a factor of 6 smaller than the spin- 
down power of the primary (assuming ju ~ 10 32 Gcm 3 one 
finds Li_2 ~ 10~ 5 L s d). This indicates that the contribution of 
this interaction to the flux of relativistic particles cannot ex- 
ceed that of the pulsar-like acceleration mechanism presented 
in this paper. Furthermore, as shown by Rafikov et al. ([1999 1, 
the energy released in the process of the interaction between 
the magnetic field and a diamagnetic body moving through the 
field is converted mainly to the thermal energy of the body and 
the energy of Alfven waves generated in this process. 

An effective acceleration of charged particles can occur 
outside the magnetosphere in a region of interaction between 
the relativistic wind ejected by the white dwarf and the mate- 
rial surrounding the system. This interaction leads to a forma- 
tion of a shock at a distance (see e.g. de Jager & Harding 1992 
and references therein) 

r s - IO'Vm V^I^fcm, (35) 

where p_24 and are the density and the thermal velocity 
of the material surrounding the system expressed in units of 
l(T 24 gcirr 3 and lO^ms" 1 , and L 34 = L sd /10 34 erg s _1 . The 
energy of particles accelerated in the shock can be as high as 

£sh,max ~ 10 13 B-6 (/ /r s ) eV, (36) 

where fi_6 is the magnetic field strength in the shock region ex- 
pressed in units of 10~ 6 (this normalization is derived by equat- 
ing the magnetic field energy density, B 2 /8n, to the energy den- 
sity of the background material, (1 /2)pv 2 ). The radiative losses 
of these particles under these conditions are dominated by the 
synchrotron emission and the energy of the corresponding pho- 
tons lies within the region of ~ 200 MeV. Furthermore, the ra- 
diation time of these electrons is close to 10 6 yr. Therefore, the 
energy radiated by the electrons on a scale of r s is small. This 
indicates that a contribution of the shock into TeV emission of 
the system is almost negligible. 

A situation in which the white dwarf could appear as an 
ejector of very high energy particles with a rate ~ 10 34 erg s _! 
has been discussed by Meintjes & de Jager (2000 1. A key as- 
sumption of their approach is that the efficiency of the drag- 
driven propeller action by the white dwarf is too low for the 
most dense blobs to be expelled from the system. These blobs, 
therefore, are able to reach the circularization radius (which for 
the parameters of AE Aqr is r c ; lc ~ 2 x 10 10 cm) and to form 
a clumpy disk surrounding the magnetosphere of the white 
dwarf. A mixing of the disk material with the magnetic field, 
which is assumed to be governed by a turbulent diffusion and 
the Kelvin-Helmholtz instability, leads to a perturbation and 
reconnection of the field lines. The energy releases in the cor- 
responding current sheets is assumed to be transferred into the 



energy of particles accelerated in these regions (for a discussion 
see also Meintjes & Venter 2005 1. 

The energy of accelerated particles and the rate of their 
ejection evaluated within the MHD-propeller model are high 
enough for the flux of TeV emission to be comparable to that 
reported by Potchefstroom and Durham groups. Furthermore, 
this model is also helpful for an interpretation of the system ra- 
dio and mid-infrared emission (Meintjes & Venter 2005}. The 
reasons for a lack of success in searching for the high-energy 
emission of AE Aqr by the Whipple group in this case might 
be a relatively low probability to detect infrequent TeV activity 
of the star (Lang et al. 1998 1, or/and the energy dependence of 
the image selection criteria, used by the Whipple group in their 
data analysis procedure, which may lead to a significant loss of 
y-ray events in case of sources with spectra harder that that of 
Crab Nebular (for a discussion see Bhat et al. H998i . 

At the same time, a conclusion that MHD-propeller model 
provides us with a comprehensive explanation of the energy 
release process in AE Aqr seems to be rather premature. The 
model is built around an assumption about a very high rate of 
plasma penetration into the magnetic field of the white dwarf. 
In particular, the diffusion coefficient is normalized by Meintjes 
& Venter (2005} to its absolute maximum value 10 I6 cm 2 s _1 , 
which for the conditions of interest is 10 orders of magnitude 
larger than the value of the Bohm diffusion coefficient. Such 
a situation is unique in astrophysical objects. In particular, it 
has never been observed in solar flares, interplanetary space 
(see e.g. Priest & Forbes 2000 1, and accretion-driven compact 
stars (Frank et al. 1985 1. Furthermore, the time of blob dif- 
fusion into the magnetic field is evaluated to ~ 30 s. This is a 
factor of 2 smaller than the free-fall time, f ff = r 3/2 /(2GM wd ) 1/2 
at a distance of closest approach of the blobs to the white 
dwarf (~ 10 1() cm). This indicates that the rate of diffusion 
evaluated within the MHD-propeller model significantly ex- 
ceeds the rate of plasma penetration into the magnetic field 
of a neutron star evaluated by Arons & Lea (119761 for the 
case of spherical accretion and by Ghosh & Lamb ( 1979 1 for 
the case of disk accretion. Finally, the diffusion time of the 
blobs into the magnetic field within this approach is signifi- 
cantly smaller than the characteristic time of drag interaction 
between the blobs and the magnetic field, which according to 
Wynn et al. J1997I is ?d rag » fff. This rises a question about 
the structure of the Ho- Doppler tomogram expected within the 
MHD-propeller model. A simulation of this tomogram within 
the MHD-propeller model may help us to answer the question 
why the contribution of the clumpy disk into the observed to- 
mogram is negligibly small. 

The above mentioned items suggest that the theoretical de- 
velopment of the MHD-propeller model remains sofar work 
in progress. Nevertheless, the basic prediction of this model, 
namely, a flux of TeV emission from AE Aqr over the level of 
10 -12 cirT 2 s _1 can be observationally tested with the MAGIC 
(Lorenz l2004i. VERITAS (Weekes et al. 120021 and HESS 
(Hinton 2004 1 experiments. The corresponding observations 
will allow us to choose between the MHD-propeller model and 
EWD-model which suggests that the TeV flux of the system is 
below the threshold of detectors used in the above mentioned 
high-energy experiments. 
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